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Abstract 

The  consistent  problem  of  the  CO2  emissions  and  the  necessity  to  find  new  energy  sources,  are  motivating  the  scientific  research  to  use  high 
efficiency  electric  energy  production’s  technologies  that  could  exploit  renewable  energy  sources  too.  The  molten  carbonate  fuel  cell  (MCFC)  due 
to  its  high  efficiencies  and  low  emissions  seems  a  valid  alternative  to  the  traditional  plant.  Moreover,  the  high  operating  temperature  and  pressure 
give  the  possibility  to  use  a  turbine  at  the  bottom  of  the  cells  to  produce  further  energy,  increasing  therefore  the  plant’s  efficiencies.  The  basic  idea 
using  this  two  kind  of  technologies  (MCFC  and  microturbine),  is  to  recover,  via  the  microturbine,  the  necessary  power  for  the  compressor,  that 
otherwise  would  remove  a  consistent  part  of  the  MCFC  power  generated.  The  purpose  of  this  work  is  to  develop  the  necessary  models  to  analyze 
different  plant  configurations.  In  particular,  it  was  studied  a  plant  composed  of  a  MCFC  500  kW  Ansaldo  at  the  top  of  a  microturbine  lOOkW 
Turbec.  To  study  this  plant  it  was  necessary  to  develop:  (i)  MCFC  mathematical  model,  that  starting  from  the  geometrical  and  thermofluidodynamic 
parameter  of  the  cell,  analyze  the  electrochemical  reaction  and  shift  reaction  that  take  part  in  it;  (ii)  plate  reformer  model,  a  particular  compact 
reformer  that  exploit  the  heat  obtained  by  a  catalytic  combustion  of  the  anode  and  part  of  cathode  exhausts  to  reform  methane  and  steam;  and  (iii) 
microturbine-compressor  model  that  describe  the  efficiency  and  pressure  ratio  of  the  two  machines  as  a  function  of  the  mass  flow  and  rotational 
regime.  The  models  developed  was  developed  in  Fortran  language  and  interfaced  in  Chemcad®  to  analyze  the  power  plant  thermodynamic  behavior. 
The  results  show  a  possible  plant  configuration  with  high  electrical  and  global  efficiency  (over  50  and  74%). 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Different  kinds  of  fuel  cells  characterized  by  high  efficiency 
and  environmental  compatibility  are  under  development  at  the 
moment.  In  particular,  the  technology  of  molten  carbonate  fuel 
cells  (MCFCs)  is  now  at  the  stage  of  scale  up  to  commercial¬ 
ization  and  many  developers  have  shown  significant  progresses. 
MCFCs  are  planar  cells  formed  by  a  matrix  (tile)  filled  with  car¬ 
bonates  and  coupled  with  two  electrodes  where  the  following 
reactions  occur: 

H2  +  CO32-  =►  H20  +  C02  +  2e"  (A)  (1) 

]-02  +  C02  +  2e“  =►  CO32-  (C)  (2) 

On  the  anodic  side,  the  water  gas  shift  reaction  occurs  too: 
H20  +  CO  &  H2  +  C02  (A)  (3) 
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The  fuel  and  the  oxidant  gases  are  fed  separately,  and  the  tile  pre¬ 
vents  gas  crossover  and  guarantees  an  adequate  ionic  conduction 
and  electronic  insulation.  To  supply  higher  power,  MCFCs  are 
used  superimposed  in  stacks,  where  cells  are  electrically  con¬ 
nected  in  series  and  separated  from  each  other  by  bipolar  plates, 
also  accomplishing  the  gas  distribution. 

High  temperature  fuel  cells,  such  as  molten  carbonate  fuel 
cells  have  the  following  advantages: 

1  No  need  for  expensive  and  easily  poisoning  catalysts  (Pt 
based). 

2  High  operating  temperature  (about  923  K)  which  enables  effi¬ 
cient  recovery  of  residual  pressure  and  waste  heat  via  gas 
turbine  bottoming  cycle  and  cogeneration. 

But  focusing  the  attention  in  the  energy  production  process 
exclusively  on  fuel  cells  is  very  restrictive  since  they  cannot 
produce  power  on  their  own  due  to  the  fact  that  hydrogen  in 
nature  is  not  present  in  enough  quantity.  Therefore,  the  fuel  must 
be  obtained  from  primary  sources,  which  at  present  time  are 
mainly  fossil  fuels  (e.g.  oil  or  natural  gas).  For  this  reason  a 


836 


F.  Orecchini  et  al.  /  Journal  of  Power  Sources  160  (2006)  835-841 


Nomenclature 

a,b 

constant  Table  5. 

Cp 

specific  heat,  J  mol-1  K-1 

d 

channel  height,  m 

E 

Nemst  potential 

Ej 

activation  energy,  J  mol-1 

F 

Faraday’s  constant,  Cmol-1 

G 

mass  flowrate  on  compressor  and 

turbine 

h 

heat  transfer  coefficient,  kW  m-2 

K-1 

AH 

formation  enthalpy,  J  mol-1 

J 

current  density,  A  m-2 

K 

constant  in  Table  1  or  abs.  rate  in 

Table  4 

n 

gas  flow  rate,  mol  m1  s-1 

nc 

rotational  regime 

ne 

electrons  transferred  in  reactions  (1)  and  (2) 

P 

pressure,  atm 

r 

reaction  rate,  molm-2  s-1 

Rq 

electrode  internal  resistance  m2 

T^tot 

electrode  total  resistance  m2 

s 

cell  component  thickness,  m 

5 

specific  gas/solid  interface  area,  m2  m  - 

T 

temperature,  K 

V 

gas  velocity,  ms-1 

V 

voltage,  V 

A,  B, 

Cir,  D ,  F  kinetic  parameters  in 

total  resistance 

expression 

Greek  letters 

p 

compressor-expansion  ratio 

Q 

catalyst  surface 

thermal  conductivity 

Tf) 

efficiency 

P 

density 

Index 

0 

ideal  condition  adimensional 

a 

anode 

c 

cathode 

cb 

catalytic  burner 

comp 

compression 

esp 

expansion 

i 

chemical  species 

j 

reaction 

n 

component 

ref 

reformer 

s 

solid 

reforming  unit  is  needed  to  convert  primary  fuels  in  hydrogen, 
which  can  later  be  fed  to  the  fuel  cells. 

Different  studies  and  hypothesis  were  proposed  to  integrate 
high  temperature  fuel  cells  and  gas  turbine  [1-3].  In  [1],  the 
predicted  performance  for  a  220  kW,  pressurized  SOFC  directly 
coupled  with  a  micro  gas  turbine  is  reported.  The  hybrid  cycle 
shows  57%  expected  electrical  efficiency.  Ref.  [2]  reports  per¬ 
formance  predicted  for  a  20  MW,  ambient  pressure,  and  internal 


reforming  MCFC-based  power  plant  fuelled  with  natural  gas.  In 
that  case,  the  gas  turbine  is  coupled  with  the  fuel  cells  via  heat 
exchangers  and  the  expected  electrical  efficiency  of  the  system 
is  about  75%.  Similar  results  are  stated  for  1  MW  size  on  [3].  On 
the  other  hand,  the  processes  reported  in  the  literature  are  very 
different  both  for  the  fuel  cell  technology  used  and  for  the  system 
configuration  adopted.  Therefore,  a  direct  comparison  cannot  be 
made  and  also  a  comparison-based  only  on  the  expected  electri¬ 
cal  efficiency  can  be  confusing.  Furthermore,  some  extrapolation 
and  simplification  of  the  models  have  been  made  in  particular 
when  MW  ranges  are  considered.  The  aim  of  this  work  is  to 
study  a  particular  power  plant  composed  of  a  MCFC  Ansaldo 
500  kW  coupled  whit  a  microturbine  TurbecTlOO  100  kW.  Sim¬ 
ilar  jobs  are  showed  in  [4].  The  principal  difference  between 
this  job  and  [4]  are  the  models  used  for  the  MCFC  and  mod¬ 
ular  integrated  reformer  (MIR).  In  fact  the  MCFC  model  used 
for  this  job  is  a  2D  model  on  the  cell’s  plane.  The  choice  of 
this  more  complicated  model  is  owing  to  the  need  to  take  under 
control  the  local  temperature  on  the  cell’s  plane  that  have  to  be 
lower  than  983  K.  This  is  the  maximum  temperature  admitted 
to  avoid  electrolyte  evaporation  and  thermal  corrosion.  Another 
difference  is  the  model  used  for  the  reformer,  in  fact,  in  this 
job  a  kinetic  approach  was  adopted  because  it  is  reputed  that 
in  this  way  it  is  possible  to  obtain  a  better  understanding  of  the 
reformer  reactions  behavior.  A  numerical  steady- state  simula¬ 
tion  of  the  process  is  set  up  and  detailed  models  of  the  compo¬ 
nents  is  developed  in  Fortran  90  and  interfaced  with  Chemcad© 
software. 

2.  Mathematical  models 

2.7.  MCFC  model 

The  fuel  cell  model  is  two-dimensional  on  the  horizontal  cell 
plane,  while  temperature  variations  along  the  vertical  coordinate 
are  neglected.  The  model  geometry  was  divided  in  three  distinct 
zones:  a  planar  solid  zone  (S)  comprehensive  of  the  two  elec¬ 
trodes,  the  bipolar  plate  and  the  electrolytic  matrix  invested  by 
the  two  crossflow  gaseous  stream  (Anodic  A  and  Cathode  C). 
The  model  is  based  on  the  following  hypotheses:  (1)  steady  state 
conditions;  (2)  uniform  voltage  on  the  fuel  cell  plane;  (3)  adia¬ 
batic  conditions;  (4)  no  radiation  heat  exchanges  between  solid 
components  and  gas  streams;  (5)  continuous  description  of  the 
gas  flow  (distributed  into  a  number  of  discrete  channels)  in  terms 
of  a  specific  rate  of  reactants  per  unit  length  of  the  fuel  cell  side; 
(6)  fully  developed  velocity  and  temperature  profiles  in  the  gas 
streams;  (7)  plug-flow  balance  equations  for  the  gas  streams;  and 
(8)  thermodynamic  equilibrium  for  shift  reaction  at  the  anode. 
Hypothesis  (2)  is  justified  by  the  high  electrical  conductivity  of 
the  current  collectors;  assumption  (4)  is  explained  by  the  con¬ 
sideration  that  (gases  being  nearly  transparent  in  the  infrared 
region)  convective  heat  transfer  is  much  greater  than  adsorp¬ 
tion  and  radiation  between  gases  and  solid.  Moreover,  being  the 
Reynolds  number  of  the  channels  less  than  100  and  the  Peclet 
number  higher  or  much  higher  than  20,  the  plug-flow  form  of 
mass  and  energy  balances  of  the  gaseous  streams  (hypothesis  7) 
is  fully  justified. 
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Table  1 

MCFC  equations 


Mass  balance 

(Eq.  (1)) 

A 

=  n  where 

(Eq.  (2)) 

C 

=  n  where 

Energy  balance 

(Eq.  (3)) 

A 

Sr^diniCpiTa  _ 
/  v  9* 

(Eq.  (4)) 
(Eq.  (5)) 


J 


ne,F 


i 

)  +  MC  -  Ta) 

k  Pi _  Ri  / 

C  Yf'GPit  =  ffl/cPiTc  +  hc{Ts  -  Tc) 


S  hc(Ts  -  Tc)  +  ha(Ts  -  Ta) 

—  (2cond  T  Q  reaz  where  Qcon  d 


=E 


Xn  8n 


82ts  ,  d2n 


8x2 


+ 


dy2 


Q 


reaz 


=  f2r,AHj 


-  VJ 


n 


7=1 


Momentum  balance 
(Eq.  (6)) 

(Eq.  (7))  ~  gy 

Electrochemical  reactions 
(Eq.  (8)) 


dPa  _  1/  Pa  Vg 


A  ^  =  Ka 


^  dpc  _  Pcvc 


d 2 
CV> 

c  d2 


(Eq.  (9)) 
(Eq.  (10)) 


S  E  =  E0  - 
S  tftnt  =  A  e 


RTy 

neF 

B/Ts 


EL'';' 

S  V  =  E  -  RtotJ 


:  +  Rq  +  DeF/Ts 


The  equations  used  are  shown  in  Table  1 .  Mass,  energy  and 
momentum  balances  of  the  gaseous  streams  are  included  in  the 
model,  together  with  the  energy  balance  of  the  solid.  In  partic¬ 
ular,  the  rate  of  each  “/’  reaction  (r7)  has  been  calculated  on 
the  basis  of  Faraday’s  law  (reactions  (1)  and  (2))  and  shift  equi¬ 
librium  condition  (reaction  (3)).  Owing  to  the  laminar  regimen 
in  the  gas  channels,  momentum  balance  is  a  linear  function  of 
gas  velocity  (Eq.  (6)  and  (7)).  The  local  kinetics  is  evaluated  by 
subtracting  the  uniform  fuel  cell  voltage  from  the  local  thermo¬ 
dynamic  voltage  (evaluated  on  the  basis  of  the  Nernst  equation 
(Eq.  (8))).  The  voltage  losses  are  related  to  the  electrical  current 
through  the  local  value  of  the  resistance  (Eq.  (10)).  The  local 
resistance  is  evaluated  using  a  semi-empirical  model  (Eq.  (9)) 


as  deducted  by  [5],  function  of  the  temperature.  The  parame¬ 
ters  involved  in  the  model  have  been  chosen  on  the  basis  of 
the  fuel  cell  geometry  ( S ,  d ,  s),  of  the  materials  involved  (A), 
of  the  flow  regime  ( h )  and  of  experimental  data  evaluated  by 
Ansaldo  (Ka  =  22,  Kc  =  46).  Finite  difference  method  has  been 
used  for  the  numerical  solution  of  the  differential  equation  sys¬ 
tems,  together  with  relaxation  method  for  the  energy  balance 
of  the  solid,  which  is  a  Fourier  problem.  The  computer  code 
has  been  developed  in  Fortran  software  and  then  interfaced  in 
Chemcad®.  The  simulation  results  are  the  maps  of  gas  and  solid 
temperatures,  electrical  current  density,  Nernst  voltage,  polar¬ 
ization,  internal  resistance,  pressure  drops  and  compositions  of 
the  gaseous  streams  on  the  fuel  cell  plane. 

2.2.  MIR  model 

The  modular  integrated  reformer  (MIR)  is  a  particular  kind  of 
plate  reformer.  Differently  to  the  conventional  tubular  reformer 
the  heating  exchange  surface  is  a  plate  so  to  make  the  reac¬ 
tor’s  dimension  more  compact.  The  particular  geometry  allow 
to  reduce  the  reactor’s  length,  so  to  limit  the  pressure  losses  too, 
giving  the  possibility  to  the  reactor  to  work  at  350  kPa,  with 
better  reforming  efficiency.  Thanks  to  the  thin  catalyst  layer, 
this  kind  of  reactor  then  allows  faster  start-up  and  faster  loads 
changes  [6] .  Its  geometry  allows  assembling  the  plates  in  stack 
so  to  make  it  still  more  compact.  A  typical  scheme  of  this  reactor 
is  shown  in  Fig.  1. 

The  reactions  that  occur  in  the  reformer  are  showed  below: 

Ref. 


CH4  +  H20  CO  +  3H2 

(4) 

CO  +  h2o  &  co2  +  h2 

(5) 

CH4  +  2H20  <*■  C02  +  4H2 

(6) 

Burn. 

h2  +  —  02  =>-  H20 

(7) 

CH4  +  202  =►  C02  +  2H20 

(8) 

co  +  Io2  =►  co2 

(9) 

m 


LLLL 


Burner  catalyst 
Reformer  catalyst 


Burner  flow 

- 


A-A  Section 


Fig.  1.  Scheme  of  a  reformer  plate. 
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Table  2 


Arrhenius  for  rate  of  reforming 


Reaction 

A  (kmol  kgcaJ  s  1 ) 

E  (kJ  mol  1 ) 

(4) 

3.72  x  1014  Pa0'5 

240 

(5) 

5.48  x  107  Pa1 

67.13 

(6) 

8.99  x  1013  Pa0  5 

243.9 

Table  3 

Arrhenius  for  absorption  rate 

Species 

A  (Pa- 1 ) 

E  (kJ  mol  1 ) 

ch4 

6.58  x  10-9 

-38.28 

CO 

8.14  x  10“10 

-70.65 

H2 

6.05  x  10“ 14 

-82.9 

h2o 

1.78  x  1010Pa 

88.68 

Table  4 

MIR  equations 


Mass  balance  on  reformer  side 


3  n 


■gf  =  Peat  fin  with  n 


=E 


vnri 


j 


n  =  jh  (  ^ch4^h2o  -  j  x  /l  with  k\  =  Aie  Eal!RT 

n  =  (pcoPh2o  ~  A^2°2 )  x  f  with  kl  =  A2e~Ea2,RT 

Qr  =  1  +  K CoPcO  +  ^H2^H2  +  KCH4PCH4  H - H>^H2(J 

Energy  equations 

=  /ZrefC^cb  -  7fef)  + 

»'  J 

(5)^MEMZcblZcb  =  ^cb(7ref  _  7cb) 


To  describe  the  reaction  rate  typical  Arrenhius  formula  were 
used.  The  activation  energy  and  the  pre-exponential  factor  to 
complete  define  the  reaction  rate  under  Ni/Al203  catalyst  are 
deducted  from  [7]  and  are  described  in  Tables  2  and  3. 

The  following  hypothesis  was  done:  (1)  steady  state  condi¬ 
tions;  (2)  complete  combustion  on  the  burner  side;  (3)  no  NOx 
produced;  (4)  no  diffusion  limit  on  the  reformer  side;  (5)  ideal 
gas  law  whit  cp  polynomial  function  of  temperature;  and  (6) 
one  dimensional  model  on  the  progress  axis  of  reaction.  The 
equations  used  to  define  the  reactor’s  behavior  are  shown  in 
Table  4. 


2.3.  TurbecTlOO 

To  define  the  turbine  and  compressor  model  it  was  necessary 
to  build  a  mathematical  model  of  the  machine  characteristic’s 
(efficiency  versus  mass  flow  and  pressure  ratio  versus  mass  flow 
as  function  of  the  two  machine  inlet  temperature,  rotational 
regime  and  mass  flow).  The  model  used  is  shown  in  Table  5. 

2.4.  Flow -sheet  and  simulation 

Fig.  2  reports  a  general  simulated  plant  flow  sheet.  Natural 
gas,  after  a  primary  treatment,  is  mixed  with  steam,  pre-heated 


Fig.  2.  Plant’s  flowsheet. 
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Table  5 

Turbec  T100  model  used 
Compressor 

ficomp  =  c\(hc)G 2  +  C2(nc)G  +  ci(h) 

>?comp  =  [1  -c4(l  ~ncf]  (|)  (^) 

where 

Cl^l(^  a(l -b/nc)+hc(hc-b)2 

C0(n  )  =  _ l«-2bn2c) _ 

C  a(\-b/hc)+hc(hc-b)2 

9  9 

. .  x  _  abhc—b  nc 

C3  a(l-b/hc)+hc{hc-b)2 

C4  =  0.4 
Turbine 


))eSp  =  [l-f3(l-/i,)2](|)  (^f) 

(3  =  0.4 


in  a  regenerative  heat  exchanger  and  fed  to  the  reformer,  where 
CO  and  H2  are  produced.  Exit  gases  passes  again  through  regen¬ 
erative  exchanger,  set  ready  to  be  processed  at  the  anodic  side  of 
the  fuel  cells  stack.  After  electrochemical  and  gas  shift  conver¬ 
sion,  residual  concentration  of  fuel  is  still  present.  In  fact  CH4 
not  converted  in  the  reformer  pass  inert  in  the  stack  and  not  all 
H2  and  CO  is  consumed.  So  fuel,  containing  effluent,  is  mixed 


Table  6 


Input  specifications 


Stacks 

4 

Cells  per  stacks 

150 

Maximum  local  T  (K) 

983 

Steam/methane  ratio 

3 

Fuel  utilization  (%) 

75 

Anhode  inlet  T  (K) 

873 

Steam  reformer  inlet  T  (K) 

443 

Operative  pressure  (kPa) 

350 

with  part  of  cathode  exhausts  and  feed  a  catalytic  burner  where 
the  heat  necessary  for  endothermic  reactions,  taking  place  in 
reforming  process,  is  produced  and  directly  exchanged.  Hot  and 
CO2  rich  gas,  exhaust  of  the  combustor’s  side,  is  mixed  with 
fresh  air  coming  from  compressor  and  fed  to  the  cathode  where 
O2  and  CO2  are  consumed  by  the  electrochemical  reaction.  Part 
of  cathode  exhaust  gas  is  sent  to  a  turbine  for  recovering  residual 
energy,  which  is  used  to  drive  air  compressor  and  to  produce  fur¬ 
ther  electrical  energy.  The  hot  effluent  finally  releases  the  heat 
necessary  to  produce  steam  for  the  reforming  reaction  and  for 
cogeneration. 

In  order  to  have  a  reference  to  compare  the  results  of  the  dif¬ 
ferent  simulations  runs,  the  set  of  boundary  conditions  reported 
in  Table  6  has  been  selected.  They  represent  proper  operating 
conditions  for  the  present  technology  of  MCFC  stacks  and  are 


-  70000  rpm 

-  65000  rpm 

-  60000  rpm 
55000  rpm 


Fig.  4.  Turbine’s  characteristics  (/3  vs.  kg  s  1 ;  77  vs.  kg  s  !). 
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kg.  s-1 

Fig.  5.  Power  produced  vs.  methane  massflowrate. 


kg.s*1 

Fig.  6.  Efficiency  vs.  methane  massflowrate. 


kg.s*1 

Fig.  7.  Maximum  local  T  vs.  methane  massflowrate. 


included  on  the  operating  conditions  forecasted  by  Ansaldo  for 
the  “Series  500”  Demo  Plant. 

A  fixed  maximum  local  temperature  of  the  solid  parts  of  the 
cells  (983  K)  is  obtained  varying  the  input  cathode  gas  temper¬ 
ature  and  flowrate.  This  constraint  is  necessary  to  reduce  the 
electrolyte  losses  and  corrosion  phenomena.  Cathode  gas  tem¬ 
perature  and  flowrate  is  controlled  by  adjusting  the  cathode  inlet 
flowrate.  This  circumstance  is  obtained  exploiting  hot  gas  sen¬ 
sible  heat,  which  raises  the  temperature  of  the  fresh  air,  after 
mixing,  varying  the  recirculation  ratio  and  the  air  inlet.  At  the 
same  time  the  air  inlet  have  to  guarantee  the  good  matching 
between  compressor  and  turbine  maintaining  high  machines’ 
efficiency.  Another  constraint  is  the  steam  to  methane  ratio:  it 
must  be  over  3  to  avoid  the  Boudouard  reaction  that  would  pro¬ 
duce  carbon  which  inhibits  catalyst  layer. 

For  a  higher  global  efficiency,  it  is  better  to  convert  the 
energy  inside  the  fuel  cells  stack,  rather  than  in  the  turbine. 
Therefore,  high  fuel  utilization  is  desired  through  the  electro¬ 
chemical  device.  On  the  other  hand,  diffusion  phenomena  have 
to  be  avoided,  so  H2  and  CO  utilization  in  the  stack  up  to  75% 
is  selected.  This  value  is  obtained  varying  current  density  and 
methane  feeding  plant.  A  443  K  temperature  was  fixed  for  the 
steam  at  mixer  with  methane  inlet,  to  assure  after  mixing  a 
408  K  flow  temperature  (in  fact  at  350  kPa  this  is  the  necessary 
temperature  to  assure  that  steam  does  not  condense).  Consid¬ 
ering  the  constraint  explained  in  Table  6  an  air  mass  flowrate 
of  0.683  kg  s-1  was  chosen  and  a  recirculation  ratio  of  65%.  In 
Figs.  3  and  4  are  showed  the  compressor’s  characteristics  and 
turbine’s  characteristics  obtained  by  the  model. 

To  choose  the  right  methane  mass  flowrate  a  sensivity  anal¬ 
ysis  was  done  varying  the  methane  feeding  plant  and  current 
density  maintaining  fuel  utilization  at  75%. 

In  Figs.  5-7  are  shown  the  results  obtained 

As  shown  by  Figs.  5  and  6  increasing  the  methane  flow  rate 
the  power  produced  and  the  plant  efficiency  are  improved  but  as 
showed  by  Fig.  7  the  maximum  local  temperature  is  increased 
to  value  not  tolerable  by  the  fuel  cell. 

To  have  a  fuel  cell  correct  operation  it  is  necessary  there¬ 
fore  renounce  to  higher  efficiency  and  choose  lower  methane 
flowrate. 

In  Table  7  are  showed  the  input  value  chosen  and  the  results 
obtained. 


Cathode  flow  inlet 


■  963-983  K 

□  943-963  K 

□  923-943  K 

□  903-923  K 


Fig.  8.  Temperature  map  on  the  cell’s  plane. 


F.  Orecchini  et  al.  /  Journal  of  Power  Sources  160  (2006)  835-841 


841 


Table  7 


Input  chosen  and  results 


Air  (kg  s- 1 ) 

0.6833 

Methane  (kgs-1) 

0.0196 

Steam  (kgs-1) 

0.0722 

Current  density  (Am-2) 

1440 

Power  produced  (kW  AC) 

516 

Electrical  efficiency  (%) 

52.7 

Thermal  power  (kW) 

212 

Global  efficiency  (%) 

74.3 

Methane  Conversion  (%) 

88.4 

Fig.  9.  Reaction  progress  in  reformer. 


Figs.  8  and  9  show  instead  the  temperature  map  on  the  cell’s 
plane  and  the  hydrogen  produced  in  the  reformer  side.  The  map 
shows  that  the  cell’s  maximum  temperature  does  not  overreach 
983  K. 

3.  Conclusion 

The  main  steps  of  this  work  are 

•  The  development  of  a  fuel  cell  numerical  model,  implemented 
in  Fortran  90,  and  interface  whit  Chemcad©. 


•  The  development  of  a  reformer/burner  unit  model  employing 
in  Fortran  90,  and  interface  whit  Chemcad©. 

•  The  development  of  the  turbine  and  compressor  characteristic 
to  completely  define  the  behavior  of  the  two  machines. 

•  Assembling  of  a  model  for  the  complete  system  using  the 
custom  and  Chemcad©library  blocks. 

•  The  study  of  the  base  case,  driven  on  preliminary  input  spec¬ 
ifications. 

The  simulation  of  the  fuel  cell  has  showed  that  global  elec¬ 
trical  efficiency  surely  can  be  kept  easily  over  50-55%,  and  a 
cogenerative  efficiency  about  75%.  This  result  is  due  to  the  uti¬ 
lization  of  the  microturbine  that  produce  the  necessary  work 
for  compressor  and  further  energy.  At  the  same  time,  the  job 
has  showed  how  the  technologic  barriers  can  reduce  anyway 
the  plant  performance.  This  job  showed  anyway  the  symbiosis 
of  the  two  machines  (MCFC,  gasturbine)  and  gets  to  first  base 
for  ulterior  studies  aimed  to  solve  the  problems  that  instead  are 
emerged  when  the  plant  works  at  partial  load. 
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